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Abstract. We present an inverse scattering approach for computing n-peakon 
solutions of the Degasperis-Procesi equation (a modification of the Camassa- 
Holm (CH) shallow water equation). The associated non- self- adjoint spectral 
problem is shown to be amenable to analysis using the isospectral deformations 
induced from the n-peakon solution, and the inverse problem is solved by a 
method generalizing the continued fraction solution of the peakon sector of the 
CH equation. 



Degasperis and Procesi showed, using the method of asymptotic integrability, that 
the PDE 

Ut - Uxxt + {b+ l)uu x = bu x u xx + uu xxx (1) 

cannot be completely integrable unless b = 2 or b = 3. The case b = 2 is the Camassa- 
Holm (CH) shallow water equation pQ, which is well known to be integrable and to 
possess multi-soliton (weak) solutions with peaks, so called multi-peakons. Degasperis, 
Holm and Hone HUE] proved that the case 6 = 3, which they called the Degasperis- 
Procesi (DP) equation, is also integrable and admits multi-peakon solutions. They 
found the two-peakon solution explicitly by direct computation. 

The purpose of this note is to briefly describe an inverse scattering procedure 
for obtaining n-peakon solutions of the DP equation. Full details will be published 
elsewhere in a longer paper [8]. Our approach is similar to that used by Beals, Sattinger 
and Szmigielski to obtain n-peakon solutions of the CH equation |SJ|S], but the present 
case does involve substantially new features; in particular, the spectral problem is of 
third order instead of second, and consequently is not self-adjoint. 

The DP equation can be written as a system for u(x, t) and m(x,t): 

nit + m x u + 3mu x = 0, (2) 
m = u - u xx . (3) 
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As shown in 3 , this is the compatibility condition for the overdetermined linear 
system 

(d x -d^ = zm^ (4) 
ipt = [2 _1 (c - d 2 x ) + u x - ud x ] if) (5) 
for a wave function ip(x,t). The constant c is arbitrary; for our purposes c = 1 is the 
appropriate choice. 

The n-peakon solution has the form 



n 



u(x,t) = J2^k{t)e^ x - x ^, m(x,t)=Yjm k (t)6(x-x k (t)), 
fc=i 



fe=i 



(6) 



where 5 is the Dirac delta distribution. This satisfies 10 by construction, while (0) 
is satisfied if and only if the functions {x k (t), m k (t)}^ =11 which describe the positions 
and heights of the peakons, evolve according to the following system of ODE (where 
sgnO = 0): 



i=l 



m k = 2 22 m kmi sgxi(x k - x t ) e" 



(7) 



The case n = 1 is trivial: mi = constant, x± = xi(0) + mi t. Also when n = 2, the 
solution can be found by straightforward integration 

In this note we will always assume that all m k > and x± < . . . < x n . To show 
that this property is preserved by the flow, assume that it holds for some value of 

t. It can be verified directly that Mi — J2 k =i mk anc ^ ^ n ~ (jXk=i TOfe ) (llfe=i(l — 

Qx k -x k+1 y}j are constants of motion. If all m k are positive, then Mi and M n are also 

positive, which implies that there is a constant mo such that < too < mj(t) < Mi 
for all t and 1 < j • < n, and that x k (t) — x k +i{t) can never become zero; hence 
xi(t) < ... < x n {t) must hold for all t. 

Since x k = ^2 i mie~^ Xk ~ Xi ^ > m k e° > mo, we see that x k — > ±oo as t — > ±oo. 
Even more is true: the peakons scatter, that is \xj — x k \ — > co as i —> ±oo 
for all j ^ k, and the particles behave asymptotically like free particles moving 
with velocities TOfc(±oo) = lim t ^± 00 m k {t) as f ^ ±oo. Since mutual distances 
between particles grow indefinitely, the asymptotic velocities are distinct, rendering 
mj(±oo) 7^ TOfc(±oo) for all j ^ k. In this sense the DP peakons belong to the 
same class of mechanical systems as the finite Toda lattice [7] and CH peakons [S] . A 
complete proof of the scattering properties will be presented elsewhere [S]. 

Now consider equation Q in the case when to is a discrete measure as in ©. 
With the t dependence suppressed, the equation reads 



ip x (x) - ipxxx{x) = z 



' 2m k 5(x - x k ) tf)(x). 



(8) 



Let a^o = ~ 00 and x n+ i = +oo. Since if) x — if) xxx = away from the support of to, the 
wave function is piecewise given by expressions of the form 

if)(x) = A k e x + B k + C k e~ x , x 6 (x k , x k+1 ) (k = 0, 1, . . . ,n). (9) 
By |JSJ, if) and ip x are continuous at each point x k , while ip xx has a jump discontinuity 
of —2zm k if)(x k ). This gives, with / denoting the 3x3 identity matrix, 



(10) 



" A k ' 




' A k -i ' 




- e -x k 




B k 


= Sk(z) 


B k -i 


S k (z)=I-zm k 


-2 


[e Xk A,e 


Cfe 




Cfe-i 




e Xk 
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Consider the particular wave function satisfying tJj(x) = e x for x < X\; that is, 
[A , B , C ] = [1, 0, 0]. For x > x n we then have ip(x) — A n (z)e x + B n {z) + C n (z)e~ x , 
where [A n , B n , C n ] = S n (z) ■ ■ ■ S2(z)Si(z) [1, 0, 0]*, so A n , B n and C n are polynomials 
in z of degree n, with coefficients depending on mi , . . . , m n and e^ 1 , . . . , e Xn . For z = 0, 
the right-hand side of JSJl is identically zero, which gives V'( a; ) = &x f° r au x ? hence 
A„(0) = 1, fl n (0) = C„(0) - 0. 

Now impose the boundary condition on the right that ij)(x) be bounded for x > x n . 
This holds iff A n = 0; in other words, the eigenvalues of the spectral problem given 
by (JSJ together with the above boundary conditions are given by the zeros of the nth 
degree polynomial A n (z). We will see below that the eigenvalues are real (in fact 
positive) and simple. Denoting them by Ai, . . . , A n , we have A n (z) = Ylk=i0- ~~ z /\k). 

As in 0Ej it will prove useful to consider an equivalent spectral problem on the 
finite interval [—1, 1]. Let y = tanh(a;/2) and define </>(?/) by ip(x) = y^2<P{v)- This 
maps 4>{x) into a piecewise quadratic function: (f)(y) = ^(^4^(1 + y) 2 + £>fc(l — y 2 ) + 
C fe (l - y) 2 ) for (vk,Vk+i) (where y k = tanh(x fc /2), y = -1, y n+1 = 1). The spectral 
problem (JHJ with boundary conditions Bq — Cq — 0, A n — is equivalent to 

-<t>vvv{y) = zg{y)<t>{y), 0(-i) = ^,(-i) = o, 0(1) = o, (11) 

where 

g(y) = ^9k5{y -yk), gk = , 1 2 , (12) 

which generalizes the string equation approach used in [2]. At each yk, the second 
derivate has a jump: 4>yy{yk+) = 4>yy{yk—) — zgk<f>(]jk)- We define a pair of Weyl 
functions W(z) and Z(z), and let bk and Ck be the residues in their partial fractions 
decompositions: 

w(z) = = - - Bn ^ = V hk 

1 ' z<t>(l) z 2zA n (z) f^z-Xk' 

(1) _ 1 B n (z) C n (z) _ ^ c k 



(13) 



y( 7 \ = = ±_ gnW ^MfJ = 



,-0(1) 2z 2zA n (z) 2zA„(z) f^z-Xk 



where we have set Ao = (so bo = 1 and co = 1/2). We will see below that the second 
Weyl function Z(z) is actually determined by the first Weyl function W(z), a fact 
that is not obvious from the definition. 

We now derive the time evolution of the scattering data {Xj, bj} defined by the first 
Weyl function W(z), when m(x,t) evolves as described by (0). Then i/;(x,t) evolves 
according to JSJl. For x < X\, equation © shows that u x = u, so in that interval 
ip(x,t) — e~ x does indeed satisfy © for all t (with our choice of c = 1). For x > x n 
we have u x = —u, which implies that ij)(x,t) = A n (z,t)e x + B n (z,t) + C n {z,t)e~ x 
satisfies (0 in that interval if and only if 

A n = 0, B n = B n /z-2A n M+, C n = —B n M + , (14) 

where M + = J2k=i m k^ Xk ■ By computing the matrix product S n ■ ■ ■ Si it is not hard 
to see that A n (z) = 1 — Miz + . . . + (—l) n M n z n , which proves that Mi and M n are 
constants of motion, as claimed above. Moreover, by analyzing how the the coefficients 
of A n depend on positions Xj , and exploiting the scattering property of the system, it 
can be seen that as t — > oo the coefficients tend to the elementary symmetric functions 
of mx(oo) < . . . < m„(oo), implying A n (z) = lim 4 ^ l30 A n (z) = Y!k =1 (l - zm k {oo)). 
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Thus the scattering property of the DP peakons manifest itself in the spectrum of the 
Dirichlet-like problem (TJ| being real and simple. 

The evolution equations l|14|) readily imply that the scattering data flows 
according to 

A fe = constant, b (t) = 1, b k (t) = b k (0)e^ Xk (k > 1). (15) 

To see how the scattering data determines c k , and hence Z(z), we proceed as 
follows. We always have Co = 1/2. Let W{z) = —B n /2zA n — J2k=i^ k /( z ~ Afe ) 
and Z(z) — C n /2zA n — X)fc=i( c fc ~ b k )/(z — Afe). From i|10|) it follows that 
B n (z) = 2zM + + 0(z 2 ), which implies W(0) = -M+. Then O gives Z(z) = 
C n /2zA n = -M + B n /2zA n = -W{Q)W{z), so c k - b k = -W(0)b k = V'* ,6/,''., A, 
for k > 1. The polynomial C n vanishes as t — » — oo (which is again seen by 
analyzing how its coefficients depend on Xj's), hence c k — b k vanishes. By l|15[) we 
have bjb k — bj(Q)b k (0) exp(l/Aj + 1/X k )t, so integration from — oo to t yields 

n , 

c k = X k b k y j . ' (fc>l). (16) 

Finally we show that the inverse spectral problem has a unique solution, 
i.e., that the j/j's and g/s are uniquely determined by the scattering data. For 
< j < n we define (1, w 2 j-i, z 2 j-i) = j^{<i>yy,ci>yA)\y=y j + and (1, w 2 j, z 2 j) = 
^(<p, 4> y , ( t ) yy)\y=y j+i-- These quantities are analogs of remainders in the theory of one 
dimensional continued fractions. Since on the interval (yj, yj+i) the solution to (|llfl 
takes the form 

4>{y) = 4>{yj+i) + 4> y {y 3 +i) {y - Vj+i) + 4>yy{yj+i-) {y - %+i) 2 / 2 , 

we obtain the following descending fractional linear transformations for the 
remainders, where lj = yj+i — yj is the length of the interval: 



w 2j 1 w 2j 

W 2 j-i — U, z 2j-l — H 

Z 2 j Z 2 j Z 2j 

W 2 j-1 1 

w2j- 2 = , Z2j-2 = 1- zgf, 



/2 



(17) 



Z 2 j~l z 2j-l 

the iteration starts at w 2n — zW(z) 1 z 2n — zZ(z) (which are known in terms 
of scattering data) and stops at w_i,,z_i. The unknown quantities {lj,gj} are 
determined in each step from the large z asymptotics of remainders known from the 
previous step as a result of the following property: all even remainders w 2 j and z 2 j 
are O(l) at z = oo, while all odd remainders w 2 j_i and z 2 j^i are O(-) there. In 
particular, denoting by the coefficient of z~ m in the expansion of a holomorphic 
function a(z) at z = oo we obtain the recovery formulas 

_ «4 ? i 

h - joT> 9j - — ~(T) ( 18 ) 

Z 2j Z 2j-1 

Mapping this back to the original variables, we have a recursive way of deriving 
formulas for the x k s and m k s in the n-peakon solution in terms of the scattering 
data. Together with lfT3|) this gives the solution {x k (t), m k (t)} of Q. In a longer 
paper jBJ, we derive closed form expressions for these quantities, and analyze the 
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dynamics in more detail. Here we merely state the results for the three-peakon case 
(n = 3): 

x 3 (t) = \n{b 1 +b 2 + b 3 ), 

(Ai-A 2 ) 2 , , (Ai — A3) 2 t 7 , (A2-A3) 2 . , 
( } = ln A 1+ A 2 ^2+ Xl+X , b lb + A2+A3 b 2 6 3 

v ' A1&1 + A 2 6 2 + A363 

(Ai-A2) 2 (Ai-A 3 ) 2 (A2-A 3 ) 2 ^ , 7 
= In (Ai + A 2 )(A 1 + A 3 )(A 2 +A3) 0l ° 2 ° 3 



m 3 (i) 



(61 + 6 2 + 6 3 ) 2 



XTT ^A 1 A 2 6 1 fe 2 + i^^Ax Wa + ^^A2A 3 6263 



A x &? + A 2 6| + A 3 6| + j^hb, + j^hb 3 + j^b 2 b 3 ' 
with similar (but more involved) expressions for roi(i) and m 2 (t). 



References 

[1] Camassa R and Holm D 1993 An integrable shallow water equation with peaked solitons Phys. 
Rev. Lett. 71 1661-4 

[2] Degasperis A and Procesi M 1999 Asymptotic integrability Symmetry and Perturbation Theory 
(Rome, 1998) ed A Degasperis and G Gaeta (River Edge, NJ: World Scientific Publishing) pp 
23-37 

[3] Degasperis A, Holm D D and Hone A N W 2002 A new integrable equation with peakon solutions 

Theoretical and Mathematical Physics 133 1463-1474 Preprint nlin. SI/0205023 
[4] Degasperis A, Holm D D and Hone A N W 2003 Integrable and non-integrable equations 

with peakons Nonlinear Physics: Theory and Experiment (Gallipoli, 2002), vol II eds M.J. 

Ablowitz, M. Boiti, F. Pempinelli and B. Prinari (Singapore: World Scientific Publishing) pp 

37-43. frepnnt |niiii. SI/ 0209008 
[5] Beals R, Sattingcr D H and Szmigielski J 1999 Multi-pcakons and a theorem of Stieltjes Inverse 

Problems 15 L1-L4 

[6] Beals R, Sattinger D H and Szmigielski J 2000 Multipeakons and the classical moment problem 

Adv. Math. 154 229-257. 
[7] Moser J 1975 Finitely many mass points on the line under the influence of an exponential potential 

Dynamical Systems, Theory and Applications ed J. Moser Lect. Notes in Phys. 38, Springer, 

Berlin, 1975 

[8] Lundmark H and Szmigielski J 2005 Degasperis-Procesi peakons and the discrete cubic string 
Int. Math. Res. Papers (imrp.hindawi.com), to appear. 



